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Abstract 

We compute a certain class of corrections to (specific) screening lengths in strongly coupled 
nonabelian plasmas using the AdS/CFT correspondence. In this holographic framework, these 
corrections arise from various higher curvature interactions modifying the leading Einstein grav- 
ity action. The changes in the screening lengths are perturbative in inverse powers of the 't 
Hooft coupling or of the number of colours, as can be made precise in the context where the dual 
gauge theory is superconformal. We also compare the results of these holographic calculations 
to lattice results for the analogous screening lengths in QCD. In particular, we apply these 
results within the program of making quantitative comparisons between the strongly coupled 
quark-gluon plasma and holographic descriptions of conformal field theory. 
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1 Introduction 

The AdS/CFT correspondence [H [21 [3l H] has proved to be a useful tool in probing strongly coupled 
physics. The quark gluon plasma [5] formed in heavy ion collisions in the RHIC and LHC may 
be strongly coupled and it is a useful exercise to see which features can be reproduced using the 
correspondence [SI El- One of the main reasons for this possible connection is the observation that a 
wide class of holographic theories [H |9] have a small ratio of shear viscosity (77) to entropy density 
(s) 0, which is what appears to be needed to explain the RHIC and LHC data pTI]. 

It is now known that while holographic theories describing isotropic plasma^ using Einstein 
gravity have r]/s = l/47r, this is no longer true when higher curvature corrections are taken into 
account [13], [HI [15]. While most of the work in the literature has focused on higher curvature 
corrections to transport coefficients, in this paper we will turn our attention to corrections to screening 
lengths originally considered in [16]. The screening lengths are defined through spatial correlators 
of Polyakov loops. When the spatial separation of the loops is large, the fall off of the connected 
contribution is exponential in separation, namely e"''''/^, where ^ is the screening length. The mass 
corresponding to the longest screening length is called the mass gap rrigap. On the holographic side, 
this corresponds to the hghtest supergravity mode [iTl [THJ [191 [2D] that is exchanged between two 
strings stretching between the boundary and the horizon. In Einstein gravity, this happens to be 
the time-time component of the metric. The mass corresponding to the exchange of the axion field 
is the Debye mass. Depending on which supergravity mode is exchanged, there is a corresponding 
screening length. Study of pure gluon theories reveal little dependence on Nc so it may be hoped 
that the differences between the QCD results and the large Nc results are going to be small [2T1 [22] . 
The holographic screening masses actually work out to be larger than the lattice results and it 
could be hoped that finite A, l/Nc corrections would lead to a better agreement. 

The most general higher curvature corrections at leading order are not known in IIB string 
theory. Only certain special classes of corrections have been worked out [251 [211 [25] • In the large A 
and large limit, it will turn out that the known curvature corrections are sufficient to compute the 
leading perturbative correction to certain screening lengths, namely the mass gap and the screening 
mass corresponding to a vector and a spin-2 exchange as explained below. In certain supersymmetric 
theories where one has AdSs x M5 where M5 is smooth compact manifold, the leading higher curvature 
correction is i?^ and gives 1/A^/^ and X^^'^/N^ corrections [T^ 1^ . In certain instances where the 

^It has become a whole new industry to see which fluid in nature has the lowest rj/s jlOj . 
^It has been recently claimed that anisotropy can lead to non- universal rj/s |12j . 



2 



dual includes fundamental fields, the leading corrections are l/N^. coming from terms [271 EB]- 
As we will argue below, knowledge of the higher curvature terms is sufficient to compute the leading 
corrections to the screening masses corresponding to the exchange of a graviton. Corrections to 
screening masses corresponding to the dilaton, axion (Debye mass), C2, B2 are all beyond the scope 
of current technology (see the discussion in [T6]). 

In a broader context, the higher derivative modifications of the leading supergravity action will 
expand the universality class of the dual CFT, by modifying the parameters in the n-point functions 
of the CFT. In particular, B? is a natural term that would appear but only in a non-supersymmetric 
context [29]. This has been studied in recent toy models [30] where the coupling of the B? terms 
were not restricted to be small. In the context of the present work, we will work in a perturbative 
framework. A further motivation in considering such corrections is to extend the phenomenology 
program initiated in [28] and studied further in [3T| [32]. 

An outline of the paper is as follows: In section |2l we review the calculation of the screening masses 
for scalar, vector and spin-two symmetry channels in the field theory. In section 3, we compute the 
leading corrections to specific screening masses. In section 4, we compare our results with QCD. 
Two appendices contain certain details about the calculations used in the paper. 

2 Screening masses 

We will be interested in the screening masses arising due to exchange of the graviton. The graviton 
fluctuations can be split into scalar, vector and spin-2 symmetry channels [33] . For convenience, the 
screening masses in the A/" = 4 plasma arising from 2-derivative gravity is reproduced [161 120] in 
table 1. Here P, C and T represent parity, charge conjugation and time reversal and in the Hilbert 
space interpretation respectively. CT corresponds to the Euclidean time reversal. 
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Table 1: Spectrum [201 IB]- We are interested in goo,goi and gij in this paper. 

Our goal is to calculate various screening masses in the deconfined phase. Algebraically it works 
out to be equivalent to working with various components of the graviton in the soliton background 
as explained in [16]. This will also facilitate a direct comparison of numerics with [19J and [20j. The 
reason for this equivalence is following: we consider modes of the form e*^*-*^*^^^ where {E,k^) and 
(t, z^) are three dimensional vectors. We get the black hole solution by a double analytic continuation 
of the coordinates t ^ it and t — )■ it, where r is the compactified spatial dimension. So the graviton 
modes in the soliton background, that are of the form e*^*, will correspond to exponentially decaying 



3 



modes in the black hole background. With k"^ = 0, E will correspond to the screening mass. Now 
we briefly discuss the properties of AdS soliton and calculate various screening masses following the 
discussion in [19]. 

The 5-dimensional AdS soliton metric is written as 
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= T?(/(^)^^' + V^.udx''dx^') + -ijT-^dr' , (1) 
with f{r)=[l-^ 

Here r^^ydx^dx'^ = —dt^ + dx^ + dy^ and r^^i^ is the 3-dimensional Minkowski metric. Here we choose 
the notation in which indices a, h etc. will be used for the 5-dimensional spacetime and yLi, v etc. for 
3-dimensional Minkowski spacetime. This geometry is constructed by double analytic continuation of 
a planar AdS black hole and the conical singularity is removed by making the r coordinate periodic 
with periodicity (3soiiton = By double analytic continuation when we get the black hole solution, 
the periodicity Psouton can be related to the temperature of the black hole T = 1/ Psouton ■ 

Now, to determine the screening masses corresponding to a graviton exchange, we solve for the 
linearized equation of motion for the graviton on the background ([T]) . We write the perturbed metric 
as 

Qab = 9ab + ehab , (2) 

where (jab denotes the background metric ([T]) and ehab is the perturbation with e ^ 1. The background 
metric is a solution of the five dimensional Einstein-Hilbert action with a negative cosmological 
constant 

1 



2£3 





12" 







(3) 



where R is Ricci scalar, ip is the Plank length and L is AdS radius. We insert ([2]) in this 5-dimensional 
effective action and collect the terms second order in perturbations; i.e., C(e^). Using this second 
order lagrangian we find the equations of motion for the perturbations. Now we make the ansatz that 
the graviton has a solution of the form hab = Hab{r)e^^'^. Here Hab{r) has only radial dependence 
and vectors z^^ and are in 3-dimensional spacetime with fc^ = — Mq; i.e., k is the 3-dimensional 
momentum vector. We choose the rest frame: fc^ = (— Mo,0,0). Now we can solve the equations of 
motion by making proper gauge choices that will also ensure that we are classifying the modes into 
scalar, vector or spin-2 excitations. 



2.1 Scalar excitations 

The scalar excitations will be related to the diagonal polarization of the graviton and also gives the 
mass gap in the dual field theory. We begin with the ansatz that the solution is of the following 
form: 

2 

Hrrir) = -^/(r)/7o(r) , 



L2 



Hrt,{r) = ik^doir). (4) 
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Here rj^^ is the metric on the Minkowski section of the ([T]). Choosing k^^ = {Mq,0), we get the 
following simplified form of the perturbations 

hrr = -^/(r)ifo(r)e-**^°*, hu = ^6o(r)e-^^^°* , 

hn = ^ao(r)e-^^oS Kr = ^/-i(r)co(r)e-^^°* , 

hrt = -tMod{r)e-'^''\ (5) 



where now i represents the spatial coordinates of the Minkowski spacetime. As explained in [19], to 
solve the equations of motion consistently, we also require 



ao(r) 



Hoir) 



(^oir) = 3^4_^4 -^o(^) and 



do{r) 



dbo{r) r^Rl dHo{r) 



2L3M2 dr {3r^-R^)L^M^ dr 

° Mr) . (6) 



{3r^ - R^oyimi 



These conditions are obtained by adopting the following strategy to solve the equations of motion: 

• First, we make a choice of the gauge by defining = Ho{r)/2. 

• By looking at the equations of motion of htr, we see that it is satisfied if we choose co(r) = 
-2RlHo{r)/{3r'-Rl). 

• Using these ao(r) and co(r) we simplify the other equations of motion. If the equation for h^x 
is subtracted from equation for hrr, we find the expression (|6]) for doi^r). 

• Now by simplifying the rest of the equations of motions, we get the following differential equation 

nr ^o)^^ + (5r ^o)^;— + (-Sr^ + i?^)2 ^o(r) - 0. (7) 

Here, we point out that to solve the equations of motions we do not need to fix bo{r). In the beginning 
we could have fixed bo{r) and then the rest of the conditions would have followed accordingly. 

The differential equation ([7]) is solved for the value of Mq such that the solution satisfies certain 
boundary conditions. The required boundary conditions are the following: first, we impose the 
condition that HQ{r) is finite at the horizon. This condition can be implemented by simply solving 
(171) near the horizon. We find that close to r = Rq, the solution should behave like 

/^oM^l-»|^(.-«o). (8) 

The second boundary condition is that Ho{r) should fall off asymptotically as This condition 

comes from the fact that metric fluctuations should fall off precisely with the rate to yield a non- 
vanishing stress-energy tensor (Tab) in the dual field theory [34J. So asymptotically we expect that 

, , L^Ci L^C2 L^Cs 
^o(r) = ^ + ^ + ^---, (9) 
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where Ci, C2, C3 etc. are dimensionless constants. 

To solve equation ([7]) numerically we express ([7]) and ([8]) in dimensionless variables u and J^o, 
where r = uRq and A^o = L^Mq/Rq. Now we solve the differential equations iteratively and find 
the numerical value of Mq such that the solution satisfies both the boundary conditions. This is 
the shooting method. We find that screening mass for the scalar symmetry channel, in units of 
temperature, is 

Mo = 2.3367rT, (10) 

which is consistent with fi9\ 120] . 



2.2 Vector excitations 

To study the vector and spin-2 excitations, we start with the following ansatz for the graviton 
polarizations 

2 

Hab = Sab^Hir) . (11) 
Here, Eab is a constant polarization tensor and satisfies 

Bra = era = = Saf^k^ Vo , (12) 

where is the momentum vector in the Minkowski spacetime. Further, to separate the vector 
excitations, we impose the condition that the non-vanishing components of the polarization tensor 
take the following form 

e^^ = = , with k.v = and v.v = 1 . (13) 

As we have chosen k^ = {M^q, 0, 0), for convenience, we can choose f = (0, 1, 0). We write H{r) = 
Hvo{t) and with these choices, the only non-zero perturbation is 

2 

Kx = ^H,o{r)e-'^'^'K (14) 

Just to clarify the notation, we point out that in M„o and H^o^r), the subscript is for vector 
excitation and the subscript 'o' is to indicate that we are working with Einstein gravity in ([3]). We 
find that the equation of motion for this perturbation is 

.;^,5^ + i^^i^„„M^0. (15) 

Now this equation of motion is solved numerically for MyQ such that the solution satisfies proper 
boundary conditions. We find that close to the horizon r = Rq, the solution behaves like 

H,o{r) = (r - Ro) . (16) 

The second boundary condition is that, similar to OH]), the asymptotic solution should be 

L^Ci L^C2 L^Cs 
Hvoir) = ^ + ^ + ^ ■■■ . (17) 

We can again express the equations (ITSl) and ffT6|) in dimensionless variables and solve the equation 
of motion numerically. We find that in units of temperature, the screening mass for vector channel 
is 

M^o = 4.3227rT , (18) 

in agreement with I2U] . 
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2.3 Massive spin-2 excitations 



The remaining graviton fluctuations, gij, describe the massive spin-2 excitations on field theory side. 
We again start with the ansatz ( fTTI) (with H{r) = Hsoir)) and the polarization tensor satisfies the 
conditions ( 1T2|) . We impose the tracelessness condition: 

V'"'e^u = 0. (19) 

Thus, now f|TT]) will describe two independent modes, one of which will be off-diagonal, i.e., e^y = 
Eyx = 1 and otherwise Sab = 0. Another mode will be the diagonal and traceless with e^x = 
—Syy = 1, and Eab = otherwise. For convenience, we will calculate the screening mass for the 
off-diagonal polarization which does not mix with other polarizations. In this case, the only non-zero 
perturbations will be 

hxy = ^H,o{r)e-'^'^'\ (20) 

and the equation of motion for it is 

rir' - Rt)^-^ + (5r^ - R'o)^^ + rL'MlH.,{r) = . (21) 
We find that close to the horizon r = Rq, the solution is given by 

Hsoir) = I - ^^{r - Ro) . (22) 

The other boundary condition is that the solution has asymptotic behavior similar to (Q. Now we 
solve the equation of motion ( ETi) numerically for Mgo and get [121 [2D] 

Mso = 3.4047rr . (23) 



3 Quantum corrections to screening 

As we have explained, we are interested in finding out the screening mass in the Yang-Mills plasma 
corresponding to a graviton exchange in the Polyakov loop correlator on the gravity side. This 
entails expanding the IIB low energy effective action to quadratic order in fluctuations. In the 10- 
dimensional action only the metric and the 5-form flux are turned on. This simplifles the problem 
enormously since we do not have to worry about terms such as RH^ where if is a combination of 
the RR and NSNS 3-form fleld strengths. One problem that we could face is that there could be 
mixing between the metric fluctuation and the scalar or the gauge fleld. Since we do not know these 
terms at higher derivative order completely, our results could be incomplete. Now we will argue that 
such mixing will not arise. Firstly, we assume that since we are in a perturbative approximation, 
there is no level crossing, z.e., the lightest mode remains the lightest mode and so on. When we 
Fourier expand, each mode comes with a factor of e~**^'* where Mj is the eigenvalue corresponding to 
that mode. Thus integration over time is only going to allow for mixing at quadratic order between 
degenerate modes. Looking at table 1, it is thus clear that the only problematic terms could arise 
due to the mixing of the metric and the dilaton-axion at quadratic order which would effect the result 
for the spin-2 exchange. Firstly note that the dilaton is constant to leading order. It gets sourced by 
the leading higher derivative correction at 0(7). Thus any fluctuation should also be at least 0(7). 
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Thus plugging this into the action would lead to contributions at least at 0(7^). Since there is an 
SL(2,Z) symmetry, we would also conclude that the axion would contribute at 0(7^). We can do 
somewhat better than this. Since there is an SL(2,Z) symmetry in IIB string theory, a term of the 
form d'^hd^(f) (where h and denote metric and dilaton fluctuations around the background) would 
necessitate the existence of d"^hd"'a where a is an axion-fiuctuation. The 10- dimensional origin of 
the latter term could be of the form CabcdC°''^'^j:C'"^^gV'^'V^C^'^\ However such terms cannot arise in 
perturbation theory due to parity conservation. As the problematic terms alluded to above will not 
arise and we need to simply focus on the curvature corrections. At this point, we must point out 
that although we are guided by type IIB string theory, the phenomenological program in [28] is more 
general and we expect corrections starting at four derivative order with unknown mixings between 
various fields. Nevertheless, if we work in a perturbative approximation, our comments above will 
also apply to the general case. 

Now the most well known curvature correction in IIB is the term with a specific contraction 
between 4 Weyl tensors dictated by supersymmetry. In general, we may expect that due to addition 
of fundamental matter, the corrections would begin at order and in a non-supersymmetric plasma 
there will also be terms and a more general set of terms. We will not consider terms involving 
covariant derivatives of the Weyl tensor in the non-supersymmetric case-this is an assumption for 
simplicity. The various tensor contractions of the Weyl tensors are shown in table 2. In the next 
section we will show in detail the analysis with which is relevant for the screening masses in the 
A/" = 4 plasma. 
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Table 2: Independent contractions of Weyl tensors in five dimensions 



Since we wish to consider general curvature corrections as in table 2, it is instructive to review the 
supersymmetric case first in detail to set out the procedure for the remaining cases. We begin 
with the 5-dimensional effective action in Einstein frame including the eight derivative correction 
term I 



(f'x\/—g 



12 a 4 



(24) 



where dimensionless coupling constant [15] 7 = C(3)/8A3/^ + -\/A/384Ai'^ ^ 1 with A being the 't 
Hooft coupling and W'^ is given by 



Wt = W 



hmnk 



(25) 



This form of the corrections appears in the supergravity action and is equal to {l/2Wf + W^), where 
is the ^^''-contraction of four Weyl tensors in Table |2j 



■^These terms originate from the dimensional reduction of the 10-dimensional action consisting of the metric and 
the five-form flux [35l [HI [26] . 

*Note that if A = 6tt, Nc = 3 corresponding to as = 0.5, then 7 w 0.003. 
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We write the soliton solution as 



ds^ = ^ {f{r)dT^ + ri^^dx'^dx") + -^rr^dr^ , (26) 



with 



/(r) = /o(r)(l + 7A(r)), 
9{r) = /o(r)(l + 7^i(r)) 

and /o(r) = (l " ^) • (27) 

We plug in this solution into the action ( 12^ and find equations of motion for /i(r) and gi{r). After 
solving these coupled differential equations we find that 

Ciri2 + 360r4i?i2 _ 285Rl^ 

/iW - ,i2(,4_^4) + ^2, (28) 

where Ci and C2 are constants. To fix these constants we impose following conditions: 

• We demand that the horizon is still at r = i?o- So gi{r) and fi{r) should be regular at r = i?0; 
i.e., the numerator of the expression for gi{r) should vanish at r = Rq. We get 

= - 75R'^ 

-15Ri(5r^ + 5r^Ri - 19R^) 
and now gi[r) = — , 



Here we notice that by choosing the given Ci, the stated condition is satisfied by both gi{r) and 
Mr). 

• To fix C2 we use the following condition: to find black hole solution of the action f l24|) we just 
need to do the double analytic continuation in the metric (j26|) : i.e., t ^ it and t — ir. In this 
solution, we can extract the background metric for the dual gauge theory by going to the asymptotic 
limit giving 

ds"^ = -f{oo)dt^ + dr^ + dx^ + dy^ . (30) 

To have the speed of light to be one in the dual gauge theory, one requires /(oo) = 1 and this fixes 
the constant C2 = 0. 

Having found the soliton solution, we remove the conical singularity by compactifing the r coor- 
dinate to be periodic and this periodicity Psoiuon will be related to the temperature of the dual field 
theory 

r = -T^ = ^(1 + 157)- (31) 

Psoliton ^ 

Now to calculate the free energy of the dual gauge theory, we use standard path integral technique 
in which we identify the Euclidean action {leucj of the bulk gravity with the ratio of the free energy 
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and temperature (w/T) of the dual field theory. To render the action finite, following [SB], we use 
background subtraction to compute the free energy: 



leading to the entropy density: 



^ = ^(1 + 157). (32) 



^ = -^ = ^^(1 + 157). (33) 



3.1 Scalar excitation 

Having studied the thermodynamic properties of this perturbed solution we now turn to computing 
the correction to the mass gap. We begin with 

9ah = (jab + (^Kb , (34) 

where e ^ 1. We insert f l5^ in the effective action f l2^ and find the terms of the lagrangian that 
are second order in e and use these terms to find equations of motion for perturbations. Now we 
make the ansatz that the solution is of the form hab = Habe~^^'^, where k and z are three-dimensional 
vectors in the Minkowski section of metric ( |26|) . Here = {t,x,y) and k is momentum vector with 
k"^ = — M^. As we want to find the correction to the mass gap, we expect that up to first order 
M = Mq + 7M1, where Mq is screening mass for the scalar excitation given by ( ITOl) and Mi is the 
first order correction. We further make the ansatz that solution is of the form: 

Hrrir) = -^f{r)Hir), 



Hr^{r) = ik^d{r). (35) 

For convenience, we choose the frame where fc^ = (M, 0, 0) so that the solution is of the form 
hab = Habe~^^^^. These expressions are similar to (jlj) but now all the terms have first order corrections 
in 7. So we write 

H{r) = Hq + 'jHi, a{r) = ao(r) + 7ai(r) , 

6(r) = bo{r) + jbi{r) , c(r) = Co(r) + 7Ci(r) , 

d{r) = do{r) + ^di{r) , M = Mq + jM^, (36) 

where Hq, ao, bo, Cq and do are zeroth order solutions. The perturbations are of the following form 

hrr = -^fir)Hir)e-^''\ hu = ^b{r)e-^^\ 

K, = ^a(r)e-^^«, Kr = ^r\r)c{r)e-^^'\ 

hrt = - iMd{r)e-'^^K (37) 
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The zeroth order terms of the equations of motion (0(7°)) give the differential equation (|7]). The 
first order terms of equations of motion will be functions of Ho{r), Hi{r), ai(r), 61 (r), Ci(r), di{r), 
Mq, Ml and their derivatives. Now we explain how we solve the first order equations of motion for 
perturbations: 

• First, we eliminate the higher derivatives of HQ{r) in the first order equations of motion (0(7)) 
by using the leading order equation of motion and its derivatives. 

• Similar to the discussion in section 2, we fix the gauge by choosing ai(r) = Hi{r)/2. 

• Now, we find the expression for Ci(r) by demanding that the equation of motion for ht^ is 
satisfied. We have given the expression for ci(r) in the appendix |Al It can be seen that the functional 
dependence of ci(r) on Hi{r) is similar to the dependence of co(r) on Ho{r) in (jG]). There are extra 
terms in expression of Ci(r) that depend on HQ{r) and these terms come from the correction jW^. 

• After using this expression for ci(r), we find that the equation for hu is also satisfied. Now we 
use the equation for the perturbation hrr and find the functional form of di{r). We give the explicit 
expressions in appendix |X1 

• After using this form of di{r), the equation for hrr is satisfied and from rest of the equations of 
motion we get the differential equation for Hi{r). Similar to the leading order solution, we fix the 
gauge by fixing either ai(r) or 61 (r) and rest of the consistency conditions follow accordingly. We 
also notice that the differential equation for Hi{r) (equation (164|1 in the appendix |A]) is similar to 
([7]), but with the leading order solution Hq^v) behaving as source. 

• Equations and the equation for Hi{r) are coupled differential equations. Now, knowing 
the value of Mq from ( ITOl) . we can solve for Mi. The boundary conditions for HQ{r) is given by ([H]) 
and (Q. Similar to the boundary condition for HQ{r), we solve the equation of motion for Hi{r) 
(equation flMj) ) near the horizon r = Rq. Hi{r) near the horizon is given by fl65|) . Similar to HqIv) 
in ([9]), the second boundary condition for the solution is that asymptotically Hi[r) falls off as 

Finally, using the shooting method we find 

M = 7rT(2.336 - 139.5147) . (38) 

A similar analysis can be done for all the perturbations given in Table [2j We have show the 
results for these terms in the Tables El IH O and El 



3.2 Vector excitation 

In this section, we calculate the correction to the screening mass for the vector channel. We start 
with a perturbed metric of the form f l34p . We insert this metric into the effective action and 
find the equations of motion for perturbations by singling out the terms O(e^). Now we make the 
ansatz that hab = HabG~^^'^i where k and z are three-dimensional vectors in the Minkowski metric 
(|26|) with k"^ = — M^. We further expect that M^ = M^o + jM^i, where M^i is first order correction 
to the screening mass. We further make the ansatz that the solution is of the following form 

2 

Hab = eab^H,{r) . (39) 

Here Eab is a constant polarization tensor and satisfies (fT^ and (fT51) . The only non-zero perturbation 
is 

2 

hr. = ^H,{r)e-''''^\ (40) 
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For this perturbation, we assume that H^{r) = H^Q^r) + ■yHyi{r), where H^q is the solution of 
differential equation f|T5|) and H^ji is the first order correction. 

To find the equation of motion for Hyi{r), we plug in (1401) in the equation of motion for the 
perturbation hrx- After inserting = M^q + 'jMyi, we expand this equation of motion in powers 
of 7. We see that the zeroth order terms in 7 will give us the equation of motion for H^oif), ^-c, 
equation ( |T5l) . We can use this zeroth order equation of motion to eliminate derivatives of HyQ{r) 
from terms first order in 7. Finally, we get equation of motion for Hyi{r) that is given in the appendix 
|A] (equation fl66|) l 

The equation of motion for Hyi{r) and Hyo{r) are coupled differential equations. Now we can find 
Mui by using the value of M^o (equation (ITSl) ) and solving this set of coupled differential equations 
such that Hyi{r) satisfies proper boundary conditions. Using ( JT6l) we can find the solution for Hyi{r) 
near the horizon and it is given by f l67|) . 

The second boundary condition for Hyi{r) is that it falls off asymptotically similar to Hyo{r), 
i.e., a function similar to f|T7|) . Now we can solve the coupled differential equations f lTSj) and f l66|) 
for Myi by converting it into dimensionless parameters (see discussion before (135]) ). Finally, we find 
that the screening mass for the vector channel is 

My = 7rT(4.322 - 398.3547) . (41) 



3.3 Massive spin-2 excitations 

We start with hab = Hah^'^^'^- We also assume that Mg = M^o + iMgi, where M^i is first order 
correction to the screening mass. Here M^q is given by ( I26p . Further, we make the ansatz that the 
solution is of the following form 

2 

Hab = Sab^Hsir) . (42) 

Here, Sab is a constant polarization tensor satisfying f ll2p and (IT^ . Now as before we focus on the 

off-diagonal mode and set e^y = Syx = 1 and Sab = otherwise. So in this case, only non-zero 
perturbation will be given by 

h.y = ^H,{r)e-^''^' . (43) 

We further assume that Hs{r) = Hso{r) + Hsi{r) and insert into the equation of motion for h^y 
Now we expand this differential equation in 7 and find that zeroth order terms give us equation of 
motion for Hso{r), that is given by fl2T|) . The first order terms will give the equation of motion for 
Hsi{r) and the full expression is given in the appendix Rl (equation fl68|) ). 

Similar to other screening masses, knowing the value of Mgo, we solve the coupled differential 
equations ( l68l) and ( 12T]) for Msi. The boundary condition for Hgo is that near the horizon, Hgo is 
given by ( l22i) and it falls off as IQ. Hsi{r) also falls of as l/r"^ asymptotically. Also, we can use ( 122|) 
and find its solution near the horizon, which is give by ( 169|) . 

Given these boundary conditions and knowing the value of M^o, we can solve the coupled differ- 
ential equations iteratively for the correction Mgi. Finally we find that the screening mass for the 
spin-2 channel is given by 

Ms = 7rr(3.404 - 167.6197) . (44) 
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4 Comparison with QCD 



In this section, we tabulate the results obtained using the procedures outlined in section 3 and ap- 
pendix B for the screening masses and r]/ s. Using these results we would like to make a phenomeno- 
logical study of certain properties of the QCD plasma. By adding higher derivative corrections we 
have introduced extra parameters in the theory. In keeping with the general strategy in [28], we 
will use input from lattice QCD for the energy density, the mass gap (and Ms) to fix these parame- 
ters. Using the fixed parameters, we will "predict" the value for rj/s and M^. Of course our results 
should be taken with a grain of salt since the corrections we have added are not the most general 
(since we have ignored covariant derivatives acting on the curvature tensor at six and eight deriva- 
tive order). However, interestingly we will find that the predictions are reasonable giving hope that 
enlarging the space of couplings may eventually lead to sensible phenomenology. From the latest 
lattice calculations, the ratio of energy density with its free field limit for 4-dimensional lattice QCD 
at temperature T ^ 2Tc is {e j Sq) lattice = 0.85 — 0.90 [23 EH]- The values of screening masses for 
2-flavour Nj = 2 QCD at 2Tc are : {M) lattice /t^T = 1.68 — 1.91; vector and spin-2 symmetry channels 
are {M.,)latt^ce/^TT = 2.76 - 3.02 and (Ms) lattice /tcT = 2.48 - 2.64 [ig EH EH]. 

For comparison with lattice results, let us consider the following action with up to eight-derivative 
corrections 




12 . 

(^L^Wabcd 



^^abcd ^ l3L'^W'''""^Wa^^Wbedf + liL^W^^ 



(45) 



where a, (3 and 7^ are dimensionless coupling constants and they are systematically suppressed by 
the powers of ip/L"^. In 7j, i runs from 1 to 5 to give the general 8-derivative terms corresponding to 
table 2. This lagrangian will give a description of supersymmetric plasma if /3 = and '~fiW^ = 'yW^, 
where is given by ( l25l) . For a non-supersymmetric plasma, /3 7^ and we also do not know the 
precise form of the 7jVr/ term in this case. 



Term 


MJtiT 


M./vrT 




4.322 + 12.965a - 114.522^2 


3.404 - 6.414a + 128.620a2 




4.322 + 45.738/3 


3.404 - 25.839/3 




4.322 + 603.16471 


3.404 + 279.74971 


wi 


4.322 + 238.11972 


3.404 - 532.63472 


wi 


4.322 - 285.40373 


3.404 - 122.10373 


wl 


4.322 - 695.93674 


3.404 - 307.49374 


wl 


4.322 - 216.64175 


3.404 - 673.80375 



Table 3: Corrections to M^,, Mg 



Implicitly in displaying these results we have assumed O^a^) ~ 0(/3), 0{q^) <^ 0{(3), 0{a^) <^ 
0(7j) and 0{(3^) ^ 0(7j). These are respected by the numerical solutions we have obtained in the 
discussion. Futhermore, we have included the a^ terms to check the consistency of our numerics. 

Let us start by assuming 'jiW^ = 'yW^; we will soon comment on the general form. Now, from 
table [6l we can see that energy density for this lagrangian is given by 



e = 



2£3 



(1 + 18a + 24a2 + 6/9 + I57) . 



(46) 
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A-TTti 1 e 

/I / / / O 


iViCLoo cicXU / /I J. 


W2 

V V 


n _ Sn/ -1- 1 1 InF'^ 


2 S36 -1- 1 7 n66n' — Q7 704-n''^ 


r V 










2 336 + 266 317^i 




(1 + 83272) 


2.336 - 532.63472 


wi 


(1 + I2O73) 


2.336 - 139.51573 


wl 


(1 + 32874) 


2.336 - 272.67374 


wl 


(1 + 68875) 


2.336 - 691.21875 



Table 4: Corrections to rj/s and mass gap 



We would like to take its ratio with the energy density of a non-supersymmetric conformal plasma 
in free field limit. To find energy density for free fields, we begin with the comparison of conformal 
anomalies of a four dimensional CFT with the one calculated using holographic techniques (28l 



^3 

a = 7i'^— and c = vr^—- (1 + 8a) . (47) 

Here a and c are central charges of CFT. Further, if we consider a free massless field theory with A^i 
vectors, A^^o scalars and N1/2 chiral fermions, we find that [2^ HP] 

124iVi + lliVi/2 + 2iVo 

a = , 48 

720 ^ ^ 

12iVi + 3iVi/2 + iVo 



120 

(50) 



15(iVo + 2Ni - 2Ni/2 



2{No + 12iVi + 3A^i/2) ■ 

Here is a constant that characterize the three-point function in CFT. This constant is found to be 
related to /3 by ^4 = 4320/3. The original calculations for this term was done in [29] in absence of the 
quadratic terms and later it was argued in [31] that even in the presence of terms, the contribution 
to ^4 is of the form 0{f3'~f, 7^). So the expression for ^4 is correct so long we are considering only 
first order corrections from and terms. We are going to normalize by the energy density for 
a collection of free bosons and fermions which is given by 

where and Nf are the number of bosonic and fermionic degrees of freedom. For a non-supersymmetric 
theory A^;, 7^ Nf where Ni, = Nq + 2Ni and Nf = 2Ni/2. Now for the comparison of energy density 
(H6|) with its free field limit, we express quantities in (ISTj) in terms of bulk gravity parameters using 

(HTD - (ED): 



^0 = ^ i^No + 2N, + -iVi/2 J = —^3 (1 + 16« + 192/3) . (52) 

Now we can use fj46|) . fl52l) and results from tables [3], HI [6] to find that up to first order in /3, 7 and 
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second order in a 



M 
Iff 

ttT 



4.32 + 12.97a - 114.52a2 ^ 45.74/5 _ 394.35^ , 



2.34 + 17.07a - 97.70a2 - 13.45/3 - 139.517 , 



3.40 - 6.41a + 128.62a2 - 25.84/3 - 167.627 , 
^(1 + 2a - Sa^ - 186/3 + I57) , 
— (1 - 8a + 112a2 - 96/3 + I2O7) . 



(53) 



(54) 



(55) 



e 



(56) 



£0 



(57) 



s 



To begin with, let us compare spectrum of type IIB theory with QCD to see what do we get. In this 
case, the supersymmetric plasma only has the correction and demanding e/eo = [0.85, 0.90] lead^ 
to 7 = [0.008,0.013]. This leads to t]/s = [0.16,0.20], M/vrT = [0.47,1.10], MjnT = [-0.94,0.81] 
and Mg/irT = [1.16, 1.91]. Thus while rj/s seems to be in the right ball park as RHIC and LHC 
data show |TT], the screening masses are underestimated. In fact for the lower limit, is also 
negative which indicates the possibility of not having a mass gap in the theory. Of course there is no 
reason for anything sensible with the leading correction which is relevant for the supersymmetric 
M = A plasma. Thus we should enlarge the space of couplings. First, we could consider adding 
fundamental matter which corresponds to turning on the term. In this case, we can use the 
lattice free energy and mass gap results to fix 7, a to get M^/nT = [1.03,2.40], rj/s G [0.123,0.154] 
and Ms/nT G [1.87,2.46]. This is already an improvement! However, note that certain values of a 
and 7, which produce < M or Mg < M, are not allowed as they violate our no level-crossing 
assumption. So let us now consider turning on and more general terms which would be 
relevant for a non-supersymmetric plasma to see what we get. 

In Table [31 m and El we have mentioned the results for individual correction terms that might arise 
in the gravitational dual of non-supersymmetric field theories. Here we observe the following pattern: 
consider just the general set of terms. If the correction to the temperature is T = 7rL^/i?o(l+Q7i) 
then the corrected mass gap is M/ttT ^ (2.34 — 9.4cj7j). Here 7^ is a dimensionless coupling 
constant and Cj depends on the different terms. Note that the above relation for the mass 
gap is approximate which we will use to get a range for the other physical quantities. Further, 
we also have the energy density e/eq = 3/4(1 + Ci'ji) and the screening mass for spin-2 channel is 
Mg/ixT (3.40 — 9.3cj7j). For the screening mass in the vector symmetry channel, we find that 
My/nT = 4.32 + miCi'^i where mi varies between —4.2 and 24.0 for different cases. The corrections 
io 7]/ s have similar qualitative behaviour and we find that the numerical coefficients of Cj7j can vary 
between 8.0 and 14.9. Now using these approximate results, we can reinstate the other corrections and 
compute the range for the coupling constants. To do that, similar to f l53|) - fl571) . we can write these 
physical quantities with contribution from all the Wi terms weighted by their respective coupling 
constants 7^. In these expressions, we can define F = Cj7j. We replace 7 — t- r/15 in fl53|) . fl55l) and 
dSSj). In dMD and ([ST]), we replace 394.357 miF and I2O7 msF. Here we have mi G [-4.2, 24.0] 
and m2 G [8.0, 14.9]. With these approximations, we have reduced the five coupling constants 7^ to 
only one, namely F. Now we can see that there are uncertainties in the screening mass for the vector 
channel because of mi, so we will compare M, Ms and e/eq with the lattice calculations. Then we 
calculate and rj/s, and compare these with lattice results. 

^This value of 7 would correspond to as = 0.14 — 0.20 for Nc — 3. 
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By using the values of {M)iattice, {Ms)iattice and {e/eo)iatuce in equations ([S3D, dSS]) and we 

get: 

a e [0.00423,0.02847], (3 G [-0.00060,0.00003], T e [0.07300,0.09179]. (58) 
Using (pD in (pD and (EZD, we find that@ 

- G [0.11494,0.19049], (59) 
s 

^ G [2.14224,4.98317]. (60) 

Thus ///s is still well withii]0 experimental limits pT]. Interestingly, the range for screening mass 
in vector symmetry channel also encompasses the lattice results (My) lattice /t^T = 2.76 — 3.02. This 
gives hope that enlarging the space of couplings along the lines of [281 El] may yield physically 
interesting results. In our numerics, we find that 0(/3) ~ 0(a;^) ^ 0(7). Curiously, the dominant 
contribution comes from the 7-dependent terms. 

5 Discussion 

In this paper we have studied corrections to certain screening masses in the Yang-Mills plasma 
at strong coupling using AdS/CFT. On the gravity side the screening masses corresponded to the 
exchange of a graviton in the correlator of Polyakov loops. We were able to extract corrections to 
the mass gap which is the lightest supergravity mode. In addition we also computed corrections to 
the screening masses arising in the vector and spin-2 channels. We expanded the space of couplings 
by adding higher derivative corrections corresponding to the addition of fundamental matter and 
breaking of supersymmetry. Using lattice input for the mass gap and Ms, we were able to get 
sensible predictions for r]/s and M„. 

Our approach may give the impression that all that one needs to do to find a realistic model 
is to add new couplings and adjust them to agree with observables available to us. It would be 
somewhat unsatisfying if that was all that could be said about this program. We would like to make 
the following interesting observation at this point. Recall that in order to get positive energy fluxes, 
we need to satisfy the following constraints for any CFT in 4 dimensions 



to 2 to 2, to 8 , , 

l-4--^U>0, 1 + t4>0, 1 + I + — 14>0, 61 

3 15 6 15 3 15 

where {c — a)/c = 8a;/(l + 8a) = ^2/6 + 4t4/45 for us. Remarkably, the ranges in (|58ll satisfy all 
these three inequalities. In other words, our findings are within the theoretical limits set forth in 
[29]. That this could happen was not at all guaranteed and we have examples of parameter spaces 
(e.g., considering the lower range for M/ttT to be 1.27) which would violate the above constraints. 

Of course our numerical analysis should be taken with a grain of salt as we did not include the 
most general corrections possible. However, it will be a useful starting point for any future work 
on improving the AdS/CFT phenomenology program initiated in [2S1 EI]- For instance, it will be 
interesting to see the effect of adding a chemical potential to the corrections [IT] or in the M = 2* 
model 1421 . 



^ Dropping the terms leads to a 10% change in the lower range or rj/s and a 2% change in the upper range of 
Mv leaving everything else the same. 



^Conservatively, for a RHIC/LHC plasma, 77/s < 0.4. 
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A Gauge conditions and other functions for eight derivative 
correction 



For the eight derivative correction ( ]25ll . the variables mentioned in (13611 are: 

- 225r^^Rf^ - 1176L^M^r^^Rl - 2808r^'^R^o + 632L^M^r^Rl^ 
+ 876r^Rl^ - lOAL^MyRl^ + 1797r^Rl^ - 1275i?f ) + 



(62) 



dr 

SR^jr^ - Rl){2L^Mlr\?>r^ - R^) - 3(36rg - 95r^i^^ + 49/?g)) 

L6rii(3r4-i?4)2 ' 

dhjr) r^Rt dH,{r) 12rX^ 

'^""^ 2L3M2 dr (3r4 - i?4)^3]Vf2 dr {ir^ - Rlfim^ '^""^ 

+ ^°'-) L^M3.n(_t.^^y (- 2iH/>'«^(3.' - Rtf 

+ 6L6Mir^''(-3r^ + R^f + L^MoV2/?^(-1944ri6 + 927r^2^^ 
+ IMlr^Rl - 843r^Rl^ + 155i?^^) - 9Mo{225r'^^ - 961r^^Rl 
+ lUr^^Rl'^ + 839r^Rlf - 80r^i2^o _ 135/^24^ 

- 32L^MoV2i?^(9r^2 + ISr^i?^ - 21r^i?^ + Si?^^) + 18Moi?o(75 r^" 

+ 551r^^R^ - mQr^^R^Q + 32Ar^Rl'^ + 7Alr^Rl^ - 315i?^°)) . (63) 
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The final equation of motion for i?i(r) is following: 

-Ro)^^ + {^r -K^)^—+ (3^4 _ ^4^2 



i^o(r) 



^6^13(3^4 _ ^4)4 

- 2219r^i?20 ^ 4057^24^ _ ^4^2^2^4^2075^24 _ 320841r20i?4 + 952785r^''i?^ 

- 1022070r^2i?J2 + 50765 7r*i?^^ - 117881r^i?f + 10291i?24)^ 

- - 2025r- - 4680LXr^^i?^ + 51111r-i?^ 
+ 6024L^M2r^°i?^ - 169776r^2^^ _ 2328L^M2r^i?J2 ^ 209412rX2 
+ 2mL^Mlr'^Rf - 94623r^i?^^ + 12573i?^°) , 

and the solution near the horizon is given by 

MoL^{48L^M^ + lOOlMoi?^ + 2MiRl){L^M^ + Rl{16L + (16 + L'^M^)Ro)) 



Hi{r) = 



+ (r - Ro) 



RliL^i + 16i?2)2 
MoL3(16 + M2L2)(48LX + lOOlMoi?^ + 2MiRl) 



ARl{Lmi + 16i?2) 
The equation of motion for Hyi (r) is following: 

+ ^^L^MiyRl + 15M^oi?o(5r^ - 91r^i?^ + 85i?^)) 

-^^(8^Xr^-36r^ + 27i.S, 

and the solution near the horizon is given by 

, , , „ , 48L«Mo^ + 2bAbL^MlRl + 2L^MoMii?g - 2880i?^ 
H^r) = -{r - Ro) LRl{Lmi + 20Rl) " 

The equation of motion for Hsi{r) is given by 



2L^Msir^^ + 48L^M>2^^ + M,oRt{75r^ - 789r^R^ + 851i?g)) 



_ ^Ed!ll^3^__3l^25r' + AQL^'MiyRl - Mt'R', - 129i?«) , 

and the solution near the horizon is 

, ^ 2(48L^M5o + 2441L^M,3oi^g + 2L^°M,^oM,ii?g + 4744M,oi^g + l&V'MsiRt) 

L^M,,Rl{Lmi + iml) 
M%{A8L^Mfo + 3977L^M2ji?2 + 2L^^MsoMsiRl + 7296i?^) 



+ (r-i?o)- 



ALmi{Lm^, + 16Rl) 
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We display below the corrections to the geometry for the various higher derivative terms considered 
in this paper 



Term 




^i(^) 




2m si^or'^m - urn) 


2i?;; 8(50r*i?^ - 87m) 




-Ur'^m + em 


2(-7r''m + ism.) 






wf 


28i?;;(5r« + br'^m - 3mr,) 


28i?^(5r« + 5r^i?^ - I9m) 


^12 


^12 


wi 


56i?^(5r« + 5r^i?^ - 3i?g) 


56i?g(5r« + 5r*i?g - 19i?«) 






wl 


15i?^(5r« + 5r^i?^ - 3i?«) 


15i?^(5r« + 5r^i?^ - 19i?g) 


^12 


^12 


wt 


29/2g(5r« + 5r^/?^ - 3i?«) 


29i?^(5r« + 5r*/2g - 19i?«) 






Wl 


74i?g(5r« + 5r^/?^ - 3/?g) 









Table 5: Corrections to metric 



B Higher derivative corrections to 77/5 

In this appendix, we review the calculation of the shear viscosity of the holographic supersymmetric 
plasma with eight-derivative correction (!25l) . We use the Kubo formula that relates the transport 
coefficients of the plasma to the field theory correlators and these correlators can be calculated using 
holographic techniques [13]. For higher derivative corrections we follow [131 

Following the prescription set in ^1], we use u = Rq/t^ which is more convenient for the hy- 
drodynamic calculations. The horizon is at m = 1 and the black hole solution for flMl) is given 
by 

ds' = -^l^dte + -^.du" + ^-{dx^ + dy^ + dr^) , (70) 

where 

f{u) = /o(w)(l + 7i/i(w)), 
9{u) = /o(u)(l +7i^i(u)) 
and /o(r) = (1 - u^) , 

fi{u) = - 15m2(-5 - + 3m^) , 

gi{u) = - 15u^(-5 - + 19m^) . (71) 

Kubo formula relates the shear viscosity to the low frequency and zero momentum limit of the 
retarded Green's function of the stress tensor 

7] = - hm -ImGg^^yiu, = 0) , (72) 
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where 

G!^y,.yi^,k = 0) = -tj dtdxdydTe'''\[T^y{x),T,y{Q)]) . (73) 

Now we translate the calculation of the correlator to the dual gravity by first calculating the effective 
action for the metric perturbation hy{t,u) = J ^|^0fc(w)e~*'^*"^*'^'^. Evaluating the action f l2^ to 
quadratic order in the fluctuation (j)k{u) yields 



1 f d^k 

^^^^ " 2£3 y J^MMu,(^)(l)l(i)~k + B{u,uj)(j)'f,(j)ii, + C{u,uj)(j)[(f)-k 

+ Diu, + Eiu)<Pl<P'U + F{u)<Pl<P'_,) + /C , (74) 

where /C is the generalized Gibbons-Hawking boundary term and its detailed expression can be found 
in [H] . Now we follow the arguments given in [H] and directly read off the shear viscosity from the 
action ^ 

V = -^Mu) + K2{u))u=l , (75) 



p 



where 



t \ r I 3uuiu) ( F'{u,uj] 
Ki{u) = lim W -— A{u, u) — B[u, u) H — 



K2(n) =lim \E{u,uj) I j ---^ 1 1 . (76) 




Note that in fl75l) we have evaluated the quantities at the horizon whereas in f l72|) it was evaluated 
at the asymptotic boundary (at u = 0). For details of the arguments about these calculations, the 
reader is suggested to refer to [41J. The functions A, B, C, D, E, F at the horizon are given by 



A{u,0) 
B{u,0) 

E{u,0) 



2R^{1 - ^2) 10R^u{l - m2)(15 + 15^2 - 29m^)7 



2\ 



3R^{1 - u , 

^ 2L5 
\mlu^{\-u^f-i 

1? ' 
16i?X(l - + 2^4)7 



F(n,0) = " ^ ^. (77) 

Using these we flnd that the shear viscosity is 

^=^^(1 + 1357). (78) 

The ratio of shear viscosity with entropy density (1331) is 

^=^(1 + 1207). (79) 
s 47r 
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Term 


\ ) X Temperature 

\ / ^ 


{ -i fvArri] X Free energy 


I i T 3TTi 1 X Tl 




/ 16 ,\ 
1 - 2« - -aM 
\ / 


(1 + 18a + 24a2) 


(1 + 10a + 136a2) 


V V 


i\ + 2/3') 




(I _ qc]R] 










wi 


(1 + 5672) 


(1 + 5672) 


(1 + 88872) 


wi 


(1 + 1573) 


(1 + 1573) 


(1 + 13573) 


wl 


(1 + 2974) 


(1 + 2974) 


(1 + 35774) 


wl 


(1 + 7475) 


(1 + 7475) 


(1 + T6275) 



Table 6: Corrections to thermodynamics and shear viscosity 



The results for other correction terms are stated in Table H] and O For the sake of completeness 
of the discussion, we have borrowed the results from [H] and [Hlj for the second order corrections to 
77. 
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